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Abstract 

A multi-configuration mixing approach built on essentially complex, symmetry- 
projected Hartree-Fock-Bogoliubov (HFB) mean fields is introduced. The mean 
fields are obtained by variation after projection. The configuration space con- 
sists out of the symmetry-projected HFB vacuum and the symmetry-projected two- 
quasiparticle excitations for even, and the symmetry-projected one-quasiparticle ex- 
citations for odd A systems. The underlying complex HFB transformations are 
assumed to be time-reversal invariant and axially symmetric. The model allows 
nuclear structure calculations in large model spaces with arbitrary two-body inter- 
actions. The approach has been applied to 20 Ne and 22 Ne. Good agreement with 
the exact shell model results and considerable improvement with respect to older 
calculations, in which only real HFB transformations were admitted, is obtained. 



1 Introduction 



The shell model configuration mixing model ([|IJ, 0, ||) in general yields a very 
good description of at least the low energy phenomena in nuclear structure physics. 
Due to the large dimensions of the configuration spaces, however, complete shell 
model calculations are restricted to rather small basis systems, typically of the size 
of the IsOd shell ||. For the description of many nuclear structure problems one 
needs much larger basis systems. Examples are the investigation of giant resonances, 
medium-heavy and heavy nuclei and even comparatively simple tasks like the study 
of negative parity states in light even-even nuclei. In all these cases one is forced 
to truncate the complete shell model expansion of the nuclear wave functions to a 
manageable number of many nucleon configurations. 

One way to achieve this is provided by the use of variational approaches, which 
leave the selection of the relevant configurations entirely to the dynamics of the 
system. The simplest models of this type are the well known Hartree-Fock (||, ||) 
and the more general Hartree-Fock-Bogoliubov (HFB) (0, ||, ||, fTOfl ) approaches. 
Though the nuclear ground state is approximated here by one single generalized 
Slater determinant only, this configuration usually accounts for a large part of the 
shell model expansion of the nuclear ground state. However, in general it breaks 
all the symmetries required by the many-nucleon Hamiltonian. Thus it cannot be 
considered as a physical state but only as some intrinsic structure, from which the 
physical components have still to be obtained with the help of projection techniques. 
Moreover, in order to obtain really optimal solutions for each set of simultaneously 
conserved quantum numbers separately, the restoration of the broken symmetries 
has to be performed before the mean field is determined by the variation. 

A whole hierarchy of such symmetry-conserving variational approaches on the basis 
of HFB-type configurations have been proposed by some of us a couple of years 
ago JO]. They have become known as the VAMPIR (Variation After Mean field 
Projection In Realistic model spaces) and the MONSTER (MOdel for handling 
many Number- and Spin-projected Two-quasiparticle Excitations with Realistic 
interactions and model spaces) approaches. 

In the VAMPIR model ([0], |13| ) the energetically lowest state of a particular 
spin-parity is approximated by a single symmetry-projected HFB vacuum and the 
underlying HFB transformation is determined by variation after the projection onto 
the desired quantum numbers. Excited states with the same quantum numbers can 
be obtained by repeating this procedure with a new HFB test vacuum which is con- 
strained to be orthogonal to all the solutions already obtained. Finally then, in this 
EXCITED VAMPIR model (QT3J, 0) the residual interaction between all the ob- 
tained solutions is diagonalized. A straightforward extension of these approaches are 
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the Few Determinant (FED) VAMPIR and the EXCITED FED VAMPIR models 
(EH) EH)> which approximate each state not by a single, but by a linear combina- 
tion of several non-orthogonal symmetry-projected HFB configurations, which are 
again determined by independent, successive variations. By such chains of varia- 
tional calculations the lowest few states of a given symmetry representation can be 
obtained, irrespective of their particular structure. 



The methods are bound to fail, however, if the complete excitation spectrum with 
respect to a particular transition operator is to be described, like e.g., in the descrip- 
tion of giant multipole resonances. If this transition operator is of one body nature, 
it is obviously preferable to consider only excited states with a similar structure like 
the corresponding ground (or yrast) state. One way to achieve this is to expand the 
nuclear wave function around a symmetry projected reference vacuum, which may 
either be the usual HFB or, e.g., a VAMPIR solution. This is the essence of the 
MONSTER approach ([11], fll7| , |4[]), in which the residual interaction is diagonal- 



ized in the space of the symmetry-projected vacuum and all the two-quasiparticle 
excitations with respect to it, if an even system is considered while for odd sys- 
tems the configuration space is limited to the symmetry-projected one-quasiparticle 
configurations. 



Unfortunately, in all applications up to now, these models had to be simplified out of 
numerical reasons. This was achieved by imposing certain symmetry restrictions on 
the underlying HFB transformations. Consequently, the corresponding HFB vacua 
do not contain all principally possible correlations, but only a particular part of 
them, which becomes more and more restricted as more symmetry requirements are 
imposed. 

In the first VAMPIR calculations only real, time-reversal invariant and axially sym- 
metric HFB transformations, which neither mix proton and neutron states nor states 
of different parity, were admitted [HJ. With this real VAMPIR approach, as such 
calculations are called in the following, only states in even-even nuclei with even 
spin and positive parity could be described. If a MONSTER calculation is based on 
such a real VAMPIR transformation, the states with different symmetries (e.g., odd 
spins) are introduced by the configuration mixing. However, odd spin, or negative 
parity states in the same even-even nucleus, or states in a neighbouring odd-odd 
nucleus may have a structure, which differs considerably from the structure of the 
real reference vacuum. Thus they cannot neccessarily be described well with the 
real MONSTER on VAMPIR approach. 



A few years ago then the VAMPIR approach has been improved ( |fL3|| , ||T8| ) by 
allowing essentially complex HFB transformations as well as parity- and proton- 
neutron- mixing. Only time-reversal and axial symmetry were kept. In this complex 
VAMPIR approach many more nucleon correlations are considered and states of 
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arbitrary spin parity in even mass nuclei can be described. 



In the present work for the first time MONSTER calculations on the basis of such 
more general complex VAMPIR transformations are performed. The corresponding 
configuration spaces are considerably larger than those of the older, more restricted 
real calculations and consequently a much better description of many states is ex- 
pected. Furthermore, now also the calculations for odd-odd nuclei can be based on 
transformations particularly derived for such systems, and only for the description 
of odd systems one has still to rely on mean fields obtained for neighbouring nuclei 
as in the older approach. However, even here the underlying transformations are 
more general and consequently the configuration spaces much larger than earlier 
and thus more correlations can be described. 

In the next section ( |2.1| ) we summarize the essential features of the VAMPIR ap- 
proach without any symmetry restrictions. We then proceed in section ( |2.2j ) by 
outlining the general formalism for MONSTER calculations on the basis of the cor- 
responding VAMPIR solutions. In section (|2.3|) then the consequences of various 
symmetry restrictions on the underlying HFB transformation are discussed and ex- 
plicit formulas for the MONSTER approach restricted to essentially complex, but 
still time-reversal invariant and axially symmetric HFB transformations are derived. 
As a first test this method is applied to the two nuclei 20 Ne and 22 Ne. Here only 
a small single particle basis, the IsOci-shell was considered. This allows to compare 
the results not only to those of the more restricted real MONSTER approach but 
also with complete shell model diagonalizations. This is done in section (|3p. Finally, 
in section the present work is summarized. 



The model space is defined by a finite, .D- dimensional set of orthonormal single 
particle states T> — \k), . . .}d- The indices i, k are standing for the set of quan- 
tum numbers characterizing the state. The corresponding creation and annihilation 
operators are denoted by {c\,c\,...}d and {q, c^, . . .}r>, respectively. They obey 
the anti-commutation relations for Fermion operators. We assume that the effective 
many-body Hamiltonian appropriate for the chosen model space is known and can 
be represented by a sum of only one- and two-body terms 
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Theory 



2.1 



The VAMPIR model 




ik ikrs 



(1) 
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where t(ik) = (i\t\k) are the one-body matrix elements of the kinetic energy (or some 
single particle energies) while v(ikrs) = (ik\v\rs — sr) denote the anti-symmetrized 
two-body matrix elements of the effective interaction. 

We then introduce quasiparticle creation- and annihilation operators and a via 
the Hartree-Fock-Bogoliubov transformation [|1(J which in matrix notation is given 
by 

^(F) \ _ p ( ct ^ _ ( A T (F) B T (F) (2) 



a(F) I c I B + (F) A+(F) \ c 



This is the most general linear transformation conserving the Fermion anti-commu- 
tation relations, provided the transformation matrix F is chosen to be unitary. Via 
the inverse transformation the Hamiltonian (|l[) can be represented in terms of the 
quasiparticle creation and annihilation operators as 

H = H°(F) + H U (F) + H 20 (F) + H 22 (F) + H 31 (F) + H A0 (F), (3) 

with the upper indices denoting the number of creators and annihilators (or vice 



versa), respectively. Explicit expressions for the various terms can be found in JTT 
The vacuum for quasiparticle annihilators is given by 

\F)=(na a (F)j |0>. (4) 

Here |0) is the particle vacuum and a enumerates the quasiparticle states. The 
product runs over all quasiparticle annihilation operators with a a (F)\0) 7^ 0. Since 
the HFB transformation generally mixes basis states of all different quantum num- 
bers, the vacuum is neither an eigenstate of the square of the angular momentum 
operator, nor of its ^-component. Furthermore it has neither good proton nor good 
neutron number. The only symmetry which is still conserved is the so called number 
parity ("HI, i.e. the vacuum contains either only even or only odd nucleon number 



components, n-quasiparticle states with respect to this vacuum \F) can then be 
defined by 

t f \F) for n = 

\F{a } n > = I ^ |F) for n = l jD (5) 



Obviously, the n-quasiparticle states, too, do violate the above mentioned sym- 
metries. In order to obtain physical states, these have to be restored before in a 
particular selection of such configurations the energetically deepest solutions are de- 
termined by variation. Since any of the n-quasiparticle configurations can be written 
as vacuum to a particular set of quasiparticle annihilators, only HFB vacua will be 
discussed in the rest of this section. 
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The restoration of the broken symmetries is achieved with the help of projection 



operators (HIT], [p~3[] ) . Good parity ir is, e.g., restored by the operator 

1 



Pi 



7T 



l + vrn . (6) 



The restoration of good proton and neutron number is equivalent to projecting the 
nuclear wave function on a good mass number A and good TVcomponent of isospin. 
The corresponding operators are is given by 

P{A) = — [d<Pe^ A S A (<P) where £ A (0) = e~^ A , (7) 
o 

and 

P(2T z ) = ^-Jd X e^S 2Tz (x) where S^{ X ) = e~ ix2f % (8) 
o 

respectively. The desired angular momentum quantum numbers can be obtained by 
an integral operator, too. It is given by 

P(IM; K) = / dQ D^ K (n)R(n) (9) 

where R(Q) is the usual rotation operator and D^j^VL) its representation in angular 
momentum eigenstates. The integration is running over the full three Euler angles. 
These projection operators commute with each other and with the Hamiltonian. In 
shorthand notation we define 

©Mif = P(IM; K)P(2T z )P(A)P(n). (10) 



Physical configurations with good symmetry S, where S represents the quantum 
numbers AT z I n , are then obtained applying the projector Q S MK on the above quasi- 
particle configurations. For any quasiparticle vacuum \F), e.g., we obtain 

\F-SM)= £ e s MK \F)f K . (11) 

K=-I 

Note that one has to take the sum over all angular momentum ^-components K. 
Otherwise the projected wave function would depend on the orientation of the in- 
trinsic reference frame 

The VAMPIR approach restricts the configuration space for the lowest (the "yrast" ) 
state of a given symmetry S to such a single symmetry projected vacuum. The 
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configuration mixing degrees of freedom /J- as well as the underlying HFB transfor- 
mation F are then determined by variation 

s _ (F-SM\H\F-SM) _ 
SE = 5 (F;SM\F;SM) ~ °" (12) 

Performing this variation one obtains the optimal description of the considered yrast 
state by only a single determinant. 

The variation leads to three sets of equations which have to be solved self-consis- 
tently. The first set, resulting from the variation with respect to the mixing degrees 
of freedom /j| in flTT|), is the diagonalization of the Hamiltonian (|I|) in the space of 
the non-orthogonal configurations Qf IK \F) 



K> 

with the constraint 



£ {H s KK ,-E s N s KK ,}f K , = (13) 



(n + N s f = i (14) 



which ensures the orthonormality of the resulting states. The Hamiltonian and the 
overlap matrices in these equations are given by 

E S KK , = (F; SMK\H\F- SMK') 

where \F; SMK) = Q S MK \F) . (15) 
N S KK , = {F-SMK\F-SMK') ' I mk\ I 



The variation with respect to the HFB transformation leads to the second set of 
equations. An elegant way to perform this variation is provided by Thouless's The- 
orem ||20|| . It states that any HFB vacuum \F d ) can be represented in terms of an 
arbitrarily chosen reference vacuum \F°), non-orthogonal to \F d ), as 




with d being an anti-symmetric DxD matrix. The quasiparticles which belong to 
the vacuum \F d ) are related to those of the vacuum \F°) via 

<F d ) ) \ -Ltd Lf ) { a(F0) ) ^ 
with 1 + d T d* = L d L+ d . 
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Consequently the variation with respect to the matrix elements of F can be replaced 
by the variation with respect to the matrix elements of d. One obtains 

ff^ = E 7 5 (L d 1 ) Q , 7 g-ys (L d x ) S p = 

where ~g j5 = J2kk> ff{F d \[H - E s ]Q s KK ,a\ (F d )a\(F d ) \F d )f s K ,. ^ 



This equation (a sort of generalized Brillouin theorem) expresses the stability of 
the solution ( |TTD against arbitrary projected two-quasiparticle states with the same 
symmetry S. 

Since the HFB vacuum is invariant under unitary transformation of the quasiparticle 
operators among themselves, these two sets of equations are not yet sufficient to 
determine the HFB transformation unambiguously As usual we use this freedom to 
diagonalize the one-quasiparticle spectrum. This yields a third set of equations 

(F\a a (F)Hal(F)\F) - 8(a, (3) (F\H\F) = E a 6(a, (3). (19) 

The eigenvalues E a are called quasiparticle energies. 



These three sets form the VAMPIR variational equations |Rj. Their solution yields 
in general already a rather good description for the yrast states. If necessary, corre- 
lating symmetry projected configurations can be obtained by successive variational 
calculations. This is done in the FED VAMPIR approach |T5|. Furthermore, using 
orthonormality contraints, the procedure can be easily extended to the description 
of excited states. This is done in the EXCITED VAMPIR [O], and the even more 
general EXCITED FED VAMPIR [[15] approaches. In principle, with the help of 
these approaches nuclear states of arbitrary complexity can be described. However, 
these methods are specifically designed, to obtain wave functions for the lowest few 
states of a given symmetry S only. They can hardly be applied if, e.g., the complete 
excitation spectrum with respect to some (usually one-body) transition operator is 
required. For such problems it is preferable to consider only specific configurations 
as they can be obtained by expanding the nuclear wave functions around a suitable 
VAMPIR vacuum. This is demonstrated in the following section. 



2.2 The MONSTER on VAMPIR approach 

For a particular symmetry S we define a configuration space 

{\Q; SMK}} = {\F; SMK}, \Faf3; SMK); a < j3; K = —I, +1} (20) 
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consisting out of the symmetry-projected vacuum (see also (|TTD) 

\F;SMK) = e s MK \F) (21) 

and the symmetry-projected two-quasiparticle states 

\Faf3; SMK) ee Q s MK ai(F)4(F)\F) (22) 

with respect to it. Note that S needs not necessarily be identical to the symmetry 
for which the underlying VAMPIR mean field is obtained. In many cases it will be 
sufficient to use just the VAMPIR transformation obtained for the ground state of 
the considered system. 

The above choice of the configuration space ensures that the total wave functions of 
the excited states being linear combinations of these (non orthogonal) configurations 



lUnSM) = \Y J \FiSMK)g s 0K]l + £ \F a(3; SMK)g^ 



1a/3K:i 

K a<f3,K 



= \Q; SMK)gQ K;i (23) 
Q 

are similar in structure as the projected vacuum and can hence easily be reached 
from the latter by, e.g., one-body transition operators. The configuration mixing 
degrees of freedom gqx-% can then be obtained by diagonalizing the Hamiltonian 
matrix ([[]) 

Yl { H QK;Q'K> ~ ESN QK;Q'K'} 9Q'K' = ( 24 ) 

Q' 

being subject to the usual ort honor mality constraint 

(gS) + N SgS = 1. (25) 

Note that the sum over Q in (|23|) and (24) includes implicitly the sum over the 
i^-components. The Hamiltonian and overlap matrix elements entering the above 
equations are given by 

H S QK ;Q>K> = (Q\ SMK\M S KK ,\Q>; SMK') 
Nqk ; q'K' = (Q;SMK\e s KK ,\Q';SMK>). 



As all shell model approaches, the above diagonalization (24 , 2q ) yields many-nucleon 
wave functions which are contaminated by spurious center-of-mass (CM) excitations 
PH] , if the single particle basis T> contains more than one major oscillator shell. In 
the MONSTER on VAMPIR approach these spurious admixtures are eliminated 
at least approximately by a method originally proposed by Giraud [EES, which is 
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based on the following argument. For a pure harmonic oscillator basis an exact 
separation of the CM motion is possible if all oscillator many-particle determinants 
with energies up to a certain nTiuj are included in the configuration space. The CM 
Hamiltonian is given in that case by 



H 



CM 



rpCM _j_ y> 



CM 



p2 l 

— h -Amuj 2 R 2 

2Am 2 



(27) 



with R = (1/A) Y.f=i r i an d P = Hf=i Pi, where r< are the spatial coordinates and p« 
the momentum operators of the particles i = 1, . . . , A. uo is the harmonic oscillator 
constant and m the nucleon mass. After diagonalizing the CM Hamiltonian in the 
nhuj space, one can identify the spurious CM components as lhu, 2hu, . . ., 
nhu excitations and can eliminate them. For any other basis it is not possible to 
choose a model space which is complete with respect to nhu oscillator excitations. 
However, the diagonalization of the CM Hamiltonian will produce energy eigenvalues 
which are still clustered around the exact excitation energies |[23| . By eliminating 
the corresponding eigenstates at least the predominant spurious components are 
removed. Therefore instead of solving equations (E3]J35|), we first diagonalize the 
CM Hamiltonian H CM in the chosen non-orthogonal basis (p0|) 



E 



CM;S 
QK;Q'K' 



rpCM;S A fS \ fS 

^ ly QK;Q'K' ) JQ'K' 







(2f 



with the constraint 

(f S ) + N s f s = 1. (29) 

The solutions with excitation energies around lhu, 2tku, . . . are considered as spu- 
rious. Out of those, denoted by |/ s ) (s — 1, . . . , n s ), we can construct a projection 
operator 

A = i-EI/«)<AI. ( 3 °) 

s=l 

which modifies the effective Hamiltonian (|I]) into 

H = P S HP S . (31) 



The diagonalization problem (24| , 25^ is then solved for the modified instead of the 
original Hamiltonian. The spurious states occur now at energies E s ~ and can 
hence be easily identified. 



2.3 Restriction to complex time-reversal invariant 
and axially symmetric HFB transformations 

If no symmetry restrictions are imposed on the HFB transformations, symmetry- 
projected vacua of the type (|TT|) can be used to describe arbitrary states in arbitrary 
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nuclei [13] . This, however, has not been achieved up to now out of numerical reasons. 
Instead, for the existing numerical realisations of the VAMPIR approaches, certain 
symmetry requirements were imposed on the underlying HFB transformations. 



So, e.g., in the first VAMPIR calculations |T2|, axial symmetry and time- reversal 



invariance were required, parity- and proton-neutron mixing were neglected and 
only real HFB transformations were admitted. Consequently such real VAMPIR 
solutions were only suitable for even spin, positive parity states in doubly even 
nuclei. Performing MONSTER calculations on top of such solutions, obviously these 
restrictions are removed and states with arbitrary spin-parity in doubly even, doubly 
odd and (in one-quasiparticle approximation) odd systems become accessible, too. 
However, for their calculation one has to rely on mean fields obtained for different 
spin values than the considered one and, for doubly odd and odd mass systems even 
for neighbouring nuclei. 



In the more recent implementations of the VAMPIR approaches then parity- as well 
as proton-neutron mixing were taken into account and essentially complex HFB 
transformations were admitted. This introduces many more correlations into the 
projected vacua as in the older calculations and furthermore makes states with ar- 
bitrary spin and parity in both doubly even and doubly odd nuclei accessible. In 
these complex VAMPIR approaches only axial symmetry and time-reversal invari- 
ance are kept. Thus the quasiparticle spectrum is still twofold degenerate and hence 



only states in even mass nuclei can be described The MONSTER approach, 
however, was up to now still limited to the use of real VAMPIR solutions. In the 
present work for the first time complex VAMPIR solutions have been used in such 
multi-configuration mixing calculations. 



The mathematical apparatus of the complex VAMPIR approach has been described 
in detail elsewhere [[n|. In the following we shall therefore only scetch the essential 
ingredients and concentrate on those features which are needed in a subsequent 
MONSTER type calculation. 

Time-reversal invariance is imposed on the HFB transformation by requiring that 
with any creator a) a also its time-reversed partner a & = fa^f -1 , where f is the time- 
reversal operator, belongs to the same quasiparticle representation. Axial symmetry 
is enforced by conserving the ^-component of the angular momentum. For the 
transformation coefficients of the HFB transformation (fj) we get in this case the 
symmetries 

A ia (F) = 5(m h m a )A ia (F) = <5(m;, m a )Af-(F), 
Bia(F) = 5(m u m a )B i5l {F) = -6(m h m a )B-* a (F), (32) 
A m (F) = and B ia (F)=0. 
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Because of these symmetries, the vacuum (|TT|) becomes time-reversal invariant and 
furthermore an eigenstate of the ^-component of total angular momentum with 
eigenvalue zero. Consequently also all the n-quasiparticle configurations (|5|) have a 
definite angular momentum z-projection. 



These assumptions simplify the projection operator (|T0| ) considerably. Now two of 
the five integrations can be performed analytically and we obtain 



O 5 



21 + 1 



d-d sin (■&)d I KK ,('&)R('&) 



(27T) 5 



dip / dx 



(33) 



with R(i9) = exp(—ii9jy) and the definitions 



w 



(^ )X )=e X p{i [|A+|(2T,)]} 



S(<p,x) = exp{- 



2^ ~ 2 



(34) 



The VAMPIR variational equations (|T2|) consisting of equations ( p!3| , pTq and [19]) are 
also simplified considerably. The diagonalization ( Jl3| ) becomes now redundant since 
the vacuum has only one fixed angular momentum ^-component K = 0. Explicit for- 



mulas for the complex VAMPIR variational equations can be found in [|T3J. Here we 
only want to mention that the complex VAMPIR vacuum contains still all possible 
two nucleon couplings. However, it does not contain all possible four- and more- 
nucleon couplings ( ||15|| , ||18|| ). This becomes clear from the basic building blocks of 
the vacuum : a natural parity four nucleon state is always built either by two 
natural parity or two unnatural parity pairs, but never by one natural parity and 
one unnatural parity pair. One natural parity pair and one unnatural parity pair 
always yield an unnatural four nucleon state. Analogously one can find the missing 
couplings for larger numbers of nucleons. States dominated by such couplings can- 
not be described well with the complex VAMPIR approach. To solve the complex 
VAMPIR variational equations one has to calculate Hamiltonian and overlap matrix 
elements in between the projected vacuum. They are also needed in the complex 
MONSTER approach and therefore shown here, although similar expressions can 
be found already in JT^]. For convenience we introduce the abbreviations 




|F A )^<j J,^ and at(F A ) = 

R(n) = R(^ip,x) = Rms(<p,x)- 



at(F) forA = l 

fla^Flfr 1 forA = 2 (35) 
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The projected overlap matrix element in between two vacua, which belong to two 
different HFB transformations F 1 and F 2 , is then given by 

7T TV 7T 

(^|0ool^ 2 ) = J ^sin(^)4(^)— f&pfdx 

^ ' -7T 

{^e [ W (^ X )(<(^^x) + vr 5 ^ 2 (^^x))]}, (36) 
where the rotated overlap has been defined as 

n 12 (n)^(F^\R(Q)\F 2 x ). (37) 

n^ 2 (f2) is calculated as described in ([0, |24j], P5|). General matrix elements in 
between two projected quasiparticle states can then be evaluated using a generalized 
version of Wick's Theorem. The only four non-vanishing elementary contractions 
are here 

(F'MF^aliF^hmFl) ee 5(a,P)ni 2 (n), (38) 
(^1^(^)^(^)1(0)1^) ee [gl 2 (Q)Un 12 (n), (39) 
(F'KiF^R^aliF 2 )^ 2 ) ee [Xr^Unl 2 ^), (40) 
(F l \kti)aUFZ)4(F 2 )\F 2 x ) ee ^(n)]^(n), (41) 



where 



Xl 2 (Q) = [Al(F 1 F 2 ;Q)]- 1 , (42) 



gl 2 (Cl) ee J B A (F 1 F 2 ;fi)Xl 2T (f2), (43) 
~g 12 (Q) = B T x (F'F 2 -n)Xl 2 (Q), (44) 

and the rotated transformation matrices are defined as 

A X (F 1 F 2 ; Q) = A + (F 1 )R(Q)A(F X ) + B + (F 1 )R*(Q)B(F X ), (45) 
B X (F 1 F 2 ; Q) = B T (F 1 )R(Cl)A(F 2 ) + A T (F 1 )R*(Cl)B(F 2 ). (46) 

Here A(F{) and B(F{) are for A = 1 just the HFB transformation matrices of the 
transformation F l ; in the case A = 2 these coefficients are multiplied with the parity 

of the basis states. R denotes the representation of R in the chosen single particle 
basis T>. Similarly the rotated energy function 

hf {n) = (F^\HR{n)\F 2 ), (47) 
which is necessary for the calculation of the projected energy matrix element 



(F'lM^F 2 ) = J d0sin(0)4(0) JL J&pjdx 

^ ' -7T 

{Ke [w&,x)(hf(#,<p,x)+7r S hf {#,<p,x))]} (4* 
12 



is obtained as 
with 

and 



~ 40* 

h\ 2 {Q) = H^F 1 ) + Hf *(F X F 2 ; Q) + 3H X (F 1 F 2 ; Q), 



H^(F'F 2 ;Cl) ee J2 H a0( Fl M 2 (^)U, 

a/3 

~ 40* 

H x (F 1 F 2 ; Q) = E^A *^ 2 ;^)]^^ 2 ^)]^, 

a/3 
7<5 



(49) 
(50) 

(51) 

(52) 
(53) 



where H°, H 20 and if 40 are Hamiltonian matrix elements of (|j) in the quasiparticle 
representation. 

In case of time-reversal invariance and axial symmetry the two-quasiparticle config- 
urations can be written as 



\FaP) = al(F)a}(F)\F), 
\FaP)=aUF)al(F)\F), 
\Faa)=al(F)aUF)\F), 



\FaB)=T\Fa(3) = ai{F)a\{F)\F), 
\F{3a) =f\FaP) = aJ(F)4(F)|F), 
with a < j3 and < m a ,m^, 



(54) 



where m a is the angular momentum ^-component of state a. Time-reversed states 
are considered explicitly. The configurations of the type \Faa) are like the vacuum 
invariant under time-reversal. The total MONSTER space contains as configurations 
\q) the set (|4|) and the vacuum 

{\q)} = {\Fap), \Faf3), \Faf3), \F0a), \Faa), \F); a<p, 0<m a ,mp}. (55) 

Out of these configurations we construct linear combinations which are either even 
or odd under time-reversal 



1 

72 



\Fa(3;SM e ) = -L B 8 MK JFaP) ± 7rs(-)'*-^e&_*JFa0) 



\FafrSM e ) ee -j= [B 8 MK JFaft ± M')'^ '^"/-/v. J'*' "0 



v af3 1 



(56) 
(57) 



where K a p = m a + mp and K a § = m a — mp. The projected states of type \Faa) 
are either even or odd depending on the spin-parity of the considered symmetry S. 
For natural spin-parity, i.e. tts{~) Is — +1> they are even, for unnatural spin-parity, 
i.e. 7r 5 (— ) Is = —1, they are odd : 



\Faa- SMe/o) 



1 +/-**(-)* ®lm\Faa) 



(51 
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The same holds for the symmetry-projected vacuum \F; SM). In shorthand notation 
we have 

The most general wave function in this space is then given by 

\MnSM) = J2{\Q;SMe)g% l + \Q;SMo}g s £} (60) 
s 

= \F;SMe/o)$'°+ £ \Faa;SMe/o)g S a $ 



(0<m a ) 



+ J2 {\ F *fr SMe)g s a ^ + \Faft SMe)g% 



a<(3 
(0<m a ,rrip) 



+ SMo)g s a ° P;t + \Fafr SMo)g%} . (61) 

The expansion coefficients are obtained by diagonalizing the Hamiltonian (or, in 
order to take care of the center of mass motion the modified Hamiltonian (|3l])) 
according to equations ( p4T , p5|) . For this purpose we order the basis states \Q) in 
such a way that first all "even" and then all "odd" configurations are listed. In this 
case the hermitean Hamiltonian matrix gets the form 



H ee H eo \ _ ( H ee i Sm H eo 



(62) 



i.e., the matrix elements between two "even" or two "odd" states become purely 
real, the mixed matrix elements purely imaginary. This is explicitly shown for the 
two-quasiparticle states in appendix (A). ([62]) can be easily brought into the form 

(63) 



1 W H ee 3m H eo \ ( 1 
-i 1 J { 3m (H eo ) T 3?e H°° ) [ % 1 



which demonstrates that only a real matrix has to be diagonalized in order to obtain 
the mixing coefficients g. For the construction of H ee , H°° and H eo it is sufficient 
to calculate the matrix elements 

(Q';SMe\H\Q; SMe), 

(Q';SMo\H\Q;SMo), (64) 
(Q';SMe\H\Q;SMo). 



Since the configurations \Q) are linear combinations of type ( p6| , p7| , p8D , this in 
turn means that we have to calculate the Hamiltonian matrix elements in between 
symmetry-projected states of type (EH). Using the explicit form of the symmetry 



14 



projector (|33|) ; these are given by 



7T 7T 



{q'\HQ s KqiKq \q) 



21 + 1 / Iff 

dtisinityd^K^j—^ j dcpj dx 



(27T) 



~^^. A )( v '|H./?(0)k/) 



+ 7rsw(^x)(g'|M(^)n|g) 
+ ™(-^>, -xWlHRiA-v, -x)\q) 



(65) 



For the overlap matrices one gets the analogous equations just replacing the Hamil- 
ton operator by the unity operator 1. Thus finally we are left with the problem to 
evaluate the rotated Hamiltonian matrix elements appearing in equation fl6"5| ) and 
the corresponding rotated overlaps. The matrix elements which have the vacuum 
on both sides are just the rotated energy ( fT7| , |l9"D and the rotated overlap (^). In 
addition there occur rotated matrix elements which have on one side the vacuum 
and on the other a two-quasiparticle state and rotated matrix elements in between 
two-quasiparticle states. A straightforward calculation yields for the former 



(F 1 |M(fi)at(F A 2 )at(F A 2 )|F A 2 ) = 

with 



Xl 2T (n)(2i/ 20 *(F 1 ) + l2H* *(F 1 F 2 -,n))xl 2 (fl) 



and 



with 



(F 1 |a /3 (F 1 )a Q (F 1 )M(fi)|F 2 ) = 

n\ 2 (Q) {h\\Q)[gl\Q)U + K(F'F^Q)U} 

[h 2 \{F 1 F 2 \ Q)] a/3 = 

2H 20 (F l ) + 2H 22 (F l F 2 ;tt) 
+ (^H U (F 1 ) - 3H 3U (F l F 2 ; n))g\ 2 (n) 
- ((H 11 ^ 1 ) -3H 31 *(F 1 F 2 ;Q))g 12 (Q)) T 
+ g 12 (Cl)(2H 20 *(F 1 ) + 12H 40 *(F 1 F 2 ;n))gl 2 (Cl) 



a/3 



(66) 



(67) 



(68) 



a/3 



(69) 
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where 



[flf(F l **;n)U = EH^siF'M^U (70) 

[Sf^^Map = Y< H ^s( Fl M 2 (m,s, (71) 

and H 11 , H 20 , H 22 , H 31 and H A0 being the matrix elements of the Hamiltonian (|3|) 
in the quasiparticle representation. The corresponding rotated overlaps are just the 
elementary contractions ( fill) and (|39|). 

The rotated Hamiltonian matrix elements in between two-quasiparticle states on 
both sides are given by 



(F 1 |a /3 (F 1 )a a (F 1 )M(fi)at(F 2 )4(F 2 )|F A 2 ) = 

[gl\n)U[g{\n)U + [xl 2 (r2)] a7 [xl 2 (fi)], 5 
-[x?(ii)U[x?(fy)H)h?((i) 



+ [X 12 (O)]^^ 1 ^ 2 ;^]^^^ 2 ^)]^^ 1 ^ 1 ^ 2 ;^)]^ 
- [Xf^U/^F^K - [X^)] Pl [h^{F'F 2 -n)U 
+ iv x (F 1 F 2 ; n)U l5 }n 12 (n), (72) 



where 



K(F^Q)} 



a/3 



+ g 12 (n) (2F 20 *(F 1 ) + 12Ff *(F X F 2 ; Clj) }xf(f2) 



a/3 



and 



MF 1 ^)], 



(73) 



J a/375 — 

E I^S^f 1 ) - E [^ 2 (^)]^[^ 2 (^)]^24^; ct (f i ) 

per ^ pi^ 

V 

[xl 2 m P ,[xl 2 (n)u. (74) 

The corresponding matrix elements of the rotated overlap are obtained as 
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(F l |o /J (F l )a (F 1 )5(n)at(F 2 )oJ(F a )|i^> = 

+ ix?(n)uix{ 2 (n)u - [x?(n)) aS [x? (n)) p ,}n?(n). (75) 

Various expressions occuring in the equations ( |66| , |68| , [T^ , and ^7], |69[) display certain 
symmetries related to time-reversal. These are presented in appendix (B) and have 
been used in the numerical calculation in order to reduce the actual number of 



matrix elements to be calculated 26 



Because of time-reversal, the vacuum as well as the two-quasiparticle excitations 
with respect to it contain only components corresponding to even mass numbers. If 
odd mass systems are to be considered, the above configuration space is therefore 
not applicable. Instead here, as already in the much more restricted real MONSTER 
calculations, we choose as configuration space the one-quasiparticle excitations of 
an even mass vacuum 

{\q)} = {\Fa),\Fa),0<m a }, (76) 
where we have introduced the definitions 

\Fa) = al(F)\F) and \Fa) = r\Fa) = 4(F) \F). (77) 

The energy spectrum can now be calculated in complete analogy to the even mass 
case. First we construct states with a good symmetry S which are either even or 
odd under time-reversal. They have the form 

\Fa; SM e ) = -L [e s Mm jFa) ± n s (-y^Q s M _ m jFcy)} , (78) 

where small letters have been used for the angular momentum jg and the angular 
momentum ^-component m a to indicate the half integer nature of these quantum 
numbers in odd mass nuclei. Denoting these configurations again by 

{\Q; SMe),\Q; SMo}} = {\Fa; SMe), \Fa; SMo}} , (79) 

the most general wave function for the odd mass case takes just the same form as 
for the even mass systems floTf ). Again the diagonalization of the Hamiltonian is 
reduced to a real problem according to equations (p2|j63| , |6^ ). 

The Hamiltonian matrix elements in between the one-quasiparticle configurations 
\q) of type (|76|) , are calculated as shown in d65|) . / as well as K q and K q i of formula 
( |65D have to be replaced by the corresponding half integers. The procedure for the 
overlap matrix elements is analogous. In the end one has to calculate the rotated 
matrix elements of in betweeen the one-quasiparticle states ([7E|). The rotated 
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Hamiltonian matrix elements are obtained as 

(F 1 |a a (F 1 )M(fi)4(F A 2 )|F A 2 ) = 

ni\Q) {hl\n)[X?(Q)U + [hf(F 1 F 2 ] Q)U} , (80) 

while the rotated overlap matrix elements are simply given by the elementary con- 
traction (0)- The symmetries of these matrix elements are displayed in appendix 
(B). In general this complex MONSTER approach for the description of odd mass 
nuclei cannot be expected to have the same quality as complex MONSTER for even 
mass nuclei. Especially when the mean field properties change quite much with the 
mass number, it is a clear drawback that the one-quasiparticle states are built on a 
complex VAMPIR determinant of a neighbouring even mass nucleus. It should be 
stressed, however, that because of the much richer structure of the complex VAM- 
PIR solutions, the resulting configuration spaces are considerably larger than in the 
real case and hence a better description than by the latter can be expected also for 
odd mass nuclei. 

Finally, we would like to mention, that for any particular symmetry S in the following 
always only one vacuum IF 1 ) = \F 2 ) = \F) will be used. Nevertheless, for the sake 
of generality, all formulas in the present section have been given for the case of two 
different vacua on both sides of any matrix element. 



3 Application to 2U Ne and 22 Ne 



As a first test the above described complex MONSTER approach was applied to 
the even-even nuclei 20 Ne and 22 Ne using only the IsOd shell as single particle basis 
T>. This rather small basis was chosen in order to enable the comparison with 
complete shell model configuration mixing calculations. The single particle energies 
t(d5/2) = -4.15 MeV, *(sl/2) = -3.28 MeV, and i(d3/2) = +0.93 MeV have 
been taken from experiment |?7| for protons as well as neutrons. As effective two- 



body interaction the mass-dependent version of the Chung and Wildenthal force 
Pj] (V(A) = V(18) x (18/A) a ) has been used, with the only difference that instead 
of a = 0.3 we took a = 1/3. This force has been adjusted by its authors to many 
experimental data in various sd nuclei and is generally accepted to be "the standard" 
force for that model space. The comparison of the complex MONSTER approach 
with the exact shell model calculation provides a severe test of its quality. In addition 
the results are compared with those of the more restricted real MONSTER approach. 

Let us first turn our attention to 20 Ne. This nucleus has two active protons and two 
active neutrons in the sd shell. The complete shell model diagonalization results 
in 640 different states : 46 I* = 0+ 97 V s = 1+ 143 r = 2+ 129 I* = 3+ 
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109 I* = 4+, 64 I* = 5+ 36 r = 6+ 12 r = 7+ and 4 J w = 8+ states. In Fig. 1 the 
spectrum is plotted up to excitation energies of 16 MeV. The complex MONSTER 
approach reproduces all the 640 states exactly. This result is even independent 
of the underlying HFB transformation, provided the latter is essentially complex 
and breaks all the symmetries except time-reversal and axiality. The configuration 
mixing degrees of freedom of the complex MONSTER approach are here sufficient 
to reproduce the complete shell model spectrum. This is understandable since the 
complex MONSTER configuration space (]59[) contains in the sd shell altogether 
277 configurations, no matter which particular (even mass) nucleus is considered. 
More precisely, one has 66 configurations of type \Faf3; SMe) and \Faf3; SMo) each, 
66 of type \F aft; SMe) and \F aft; SMo) each, 12 of type \Faa; SMe/o), and the 
projected vacuum \F; SMe/o). Out of these, 57 configurations can contribute to 
spin 1 = states, 151 are available for 1 = 1 states, 223 for 1 = 2 states, 259 for 
1 = 3 states, 275 for 1 = 4 states, and finally all 277 are available for each of the spin 
I = 5, 6, 7, . . . , 14 states. Since the shell model spectrum is exactly reproduced, we 
can conclude that for 20 Ne this configuration space does contain all the shell model 
configurations. 

The complex MONSTER approach is able to reproduce even those states, which 
cannot be described with the complex VAMPIR alone. These are states which 
contain the "missing couplings", mentioned in section (|2.3|) , as irreducible substruc- 
tures. In the case of 20 Ne, e.g., the lowest 3 + and 5 + shell model states belong 
to this category. Unfortunately there is no simple "geometrical" explanation for 
the "missing couplings". To identify them, one has to write each state in terms of 
two-particle couplings using the coefficients of fractional parentage. For bigger basis 
systems this becomes very difficult or even impossible. For four nucleons in the sd 
shell, however, this can be done easily : out of all 29 (<i5/2) 4 configurations just two 
are inaccessible by the complex VAMPIR : the [J 71- = 3 + ]t=o and the [I n = 5 + ]t=o 
with maximal seniority, both being the dominant structures in the corresponding 
yrast states (see also [0). These states are exactly reproduced by the complex 
MONSTER calculation. This is an example which shows that couplings, which are 
missing in the complex VAMPIR approach, can be accounted for by the configura- 
tion mixing of a subsequent complex MONSTER calculation. However, we would 
like to stress "missing couplings" do only occur because of the restriction to an axi- 
ally and time-reversal symmetric HFB transformation. Using a HFB transformation 
without symmetry restrictions this problem is removed in the VAMPIR approach, 
too. 



For comparison Fig. 1 also displays the spectrum obtained by the real MONSTER 
approach built on the mean field resulting from a real VAMPIR calculation for the 
+ ground state : F(0 + ). In the real case the MONSTER approach has much less 
degrees of freedom than in the complex case, and the MONSTER spectrum depends 



on the choice of the VAMPIR determinant 121. First, a real VAMPIR calculation 
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in sd shell accounts only for 20 variational degrees of freedom, corresponding to the 
number of independent HFB transformation coefficients, while the complex VAM- 
PIR has 56, i.e. almost a factor of three more degrees of freedom. Second, in a real 
MONSTER calculation there are altogether only 73 configurations available : the 
vacuum \F; SM), 12 configurations of type \Faa; SM), and 30 of type \Fa(3; SM) 
and \Fa/3; SM) each. 21 of these configurations can contribute to I = states, 30 to 
1 = 1 states, 61 to / = 2 states, 56 to / = 3 states, always 60 to the / = 5, 7, 9, 11, 13 
states, and always 73 to the I = 4,6, 8, 10, 12, 14 states. These are only about one 
fourth of those available in the complex MONSTER approach. Therefore it is to 
be expected that the energies in 20 Ne obtained with the real MONSTER approach 
are higher than the exact ones, which are reproduced by the complex MONSTER 
calculation. It can be seen, however, that the difference in energy is small. On 
average one has for the yrast states a deviation of ~ 374 keV, i.e. ~ 0.9 % compared 
to the total energy of the shell model ground state, which is —41.412 MeV. The 
largest deviation, 624 keV, occurs for the I + = 5 + state, the smallest for I + = 6 + 
and amounts to 135 keV. 

Finally we want to emphasize that the exact reproduction of the shell model spec- 
trum of 20 Ne by the complex MONSTER approach is a very stringent non-trivial 
test of the computer code. 

As second example 22 Ne was studied. For this nucleus there are in total 4206 shell 
model configurations : 216 r = 0+, 534 r = 1+, 777 V = 2+, 798 r = 3+ 
723 r = 4+, 525 r = 5+ 345 r = 6+, 177 r = 7+, 81 r = 8+, 24 r = 9+, and 
6 F = 10 + ones. In Fig. 2 for each spin only the lowest states resulting from the 
shell model diagonalization are plotted. The energies obtained by the complex and 
the real MONSTER approaches, respectively, are plotted up to 16 MeV excitation 
energy. For both MONSTER calculations the underlying HFB transformation was 
determined by a corresponding VAMPIR calculation for the + ground state. On 
average the complex MONSTER reproduces the shell model energies of 22 Ne rather 
well, but for most spin values not exactly, as it was the case in 20 Ne. Here, for the 
lower spin values I n = + , . . . , 7 + not all the above listed shell model configurations 
are contained in the MONSTER configuration space and only the states with angular 
momenta I w = 8 + , . . . , 10 + are reproduced exactly. Small deviations are even found 
for the 7 + states though the number of quasiparticle configurations is here larger than 
the correponding shell model dimension. This is due to the non orthogonality of the 
symmetry-projected determinants which in general causes some linear dependencies 
inbetween the various MONSTER configurations. 

Again the complex MONSTER approach yields considerable improvements with 
respect to its real counterpart : the yrast states are more bound and many more of 
the excited shell model states can be described than in the latter case. 
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In Fig. 3 the energies of the yrast states in 22 Ne as obtained by various approaches are 
compared. The leftmost column shows the exact energies obtained by the complete 
shell model diagonalization. Then, from left to right, the energies calculated with 
the complex MONSTER, the complex VAMPIR, the real MONSTER, and the real 
VAMPIR are presented. Both MONSTER calculations have been based on the 
corresponding VAMPIR transformations obtained for the + ground state. Note 
that with the real VAMPIR approach only even spin states can be described. The 
average deviation of the energies obtained by the complex MONSTER approach from 
the shell model results is 103 keV. That is about 0.2 % of the shell model ground 
state energy of -58.820 MeV. The yrast states resulting from the real MONSTER 
approach display an average deviation of 539 keV from the exact results. This 
amounts to about 0.9 % of the shell model ground state energy. 

Except for the + ground state, where, because of the stability of the complex 
VAMPIR solution with respect to arbitrary symmetry-projected two-quasiparticle 
excitations with K = 0, both methods yield identical results, and except for the 9 + 
and 10 + state, where the latter approach already reproduces the exact energies, the 
complex MONSTER solutions are more bound than the corresponding VAMPIR 
ones. This is essentially due to the admixture of configurations with Fur- 
thermore it becomes obvious that the odd spin yrast states 1 + , 3 + , 5 + and 7 + are 
described rather poorly by complex VAMPIR alone. One can conclude that these 
states are again dominated by components corresponding to the above discussed 
"missing couplings". In fact, whenever the MONSTER result for an yrast state 
is considerably more bound than the correponding VAMPIR result this conclusion 
holds. This opens the possibility to identify "missing couplings" even in large model 
spaces where a decomposition of the configurations with the help of the coefficients 
of fractional parentage is not possible. 

The even spin states obtained with the real VAMPIR approach are on average 
360 keV less bound than the corresponding complex VAMPIR solutions. With the 
real MONSTER approach this average deviation is reduced to 98 keV. Exceptions 
are again the ground state, where MONSTER and VAMPIR solutions are again 
identical, and the 10 + state, which is exactly reproduced by the real VAMPIR, 
while the MONSTER calculation yields 101 keV less energy. This indicates that 
configuration mixing on top of the real mean field determined for the + ground 
state is not suitable to describe this high spin excitation. 

Finally we investigated the dependence of the complex MONSTER energies on the 
mean field determined by the preceeding complex VAMPIR calculation. For this 
purpose we have built always the full complex MONSTER spectrum on each of the 
yrast solutions obtained with the complex VAMPIR approach. The same was done 
in the real approximation. Obviously in this case only mean fields obtained for even 
spin states could be used. 
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As a typical example for the results of this investigation, Fig. 4 displays the de- 
pendence of the energies of the five lowest + states in 22 Ne on the spin-value used 
in the preceeding VAMPIR calculation. In addition, the leftmost column presents 
the energies of the three lowest shell model states. Solid lines denote the results of 
the complex MONSTER, dotted ones those of the real MONSTER approach. As 
can be seen, the real MONSTER results do depend only weekly on the underlying 
transformation. Here the differences in the various mean fields seem to be not very 
specific and can hence be compensated almost entirely by the configuration mixing. 
For the complex MONSTER results this does only hold as long as the + , 2 + or 4 + 
VAMPIR transformations are taken. The higher even spin transformations, how- 
ever, become now inadequate for the description of the + states. Thus one may 
conclude that the complex VAMPIR transformations display a much stronger de- 
pendence on the considered angular momentum than the more restricted real ones. 
As discussed above, due to "missing couplings", the odd spin yrast states are only 
poorly described by the complex VAMPIR approach. This is also reflected in the 
bad description of the + states, if the MONSTER diagonalization is based on an 
odd-spin transformation. 

Similar dependencies on the underlying transformations as for the + states are 
obtained for other low spin values, too. Only for very high angular momenta where 
the complex MONSTER approach exhausts the complete shell model spaces the 
results become independent of the underlying transformation. 

Finally, we would like to stress that in a MONSTER calculation (except for spin + ) 
the lowest yrast solution for a given angular momentum is not necessarily obtained 
using the VAMPIR transformation for this particular spin, but may result from a 
mean field derived for a neighbouring spin value. This was already found in real 



MONSTER calculations [12] and is due to the fact that mean field and configuration- 



mixing degrees of freedom are varied successively and not simultaneously. 



4 Summary 



In the last decade a couple of variational methods have been developed which have 
become known as the VAMPIR family. They all work with symmetry-projected HFB 
quasiparticle vacua as test wave functions, differ, however, in the degree of sophis- 
tication by which the underlying HFB transformations as well as the configuration 
mixing is determined. Originally, out of numerical reasons, in these methods only 
rather restricted HFB transformations were admitted : time-reversal invariance and 
axial symmetry were required, parity- and proton-neutron mixing were neglected 
and only real transformation coefficients were used. In the meantime most of these 
restrictions have been removed. Only the requirement of time-reversal and axiality 
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was kept and, e.g., essentially complex transformations were allowed. By these gen- 
eralisation many more correlations are accounted for and the range of applicability 
of these few determinant approaches was considerably increased. 

Unfortunately, by construction these methods are restricted to the lowest few states 
of a given symmetry only. Very often, however, one is interested in a complete set 
of particular excitations, e.g., those with respect to a specific transition operator. 
For such problems it is preferable to use a multi-configuration instead of a few de- 
terminant description. This is done in the MONSTER approach, which expands the 
nuclear wave function around the VAMPIR solution for the ground (or any other 
yrast) state and obtains the excitation spectrum by diagonalizing the chosen Hamil- 
tonian in the space of the latter and all the symmetry-projected two-quasiparticle 
configurations with respect to it. Up to now such MONSTER calculations have 
only been done using the severly symmetry-restricted real HFB transformations 
mentioned above. In the present work, now the mathematical formalism as well as 
the numerical realisation of the MONSTER approach on the basis of the much more 
general complex HFB transformations has been developed. 

With respect to the previous real implementation the complex MONSTER approach 
has a couple of essential advantages : first of all, the underlying HFB transformations 
account already for many more correlations than those used in the earlier calcula- 
tions. Second, for the description of doubly odd systems one has not any more to 
rely on mean fields obtained for neighbouring doubly even nuclei but can use HFB 
transformations obtained for the particular system under consideration. Only for 
odd mass systems, which because of the even number parity of the underlying ref- 
erence vacuum have still to be described by symmetry-projected one-quasiparticle 
states, one needs transformations from neighbouring nuclei. Third, the configura- 
tion spaces become considerably larger than in the more restricted real case and 
consequently a more complete and more detailed description of the nuclear spectra 
is to be expected. Last but not least, the admixture of two-quasiparticle configura- 
tions enables the description even of such states which are dominated by structures 
which are not contained in the symmetry-projected, time-reversal invariant HFB 
vacua used as test wave functions in the complex VAMPIR approach. 

As a first test the complex MONSTER approach was applied to the nuclei 20 Ne and 
22 Ne. For these calculations the IsOd shell was chosen as basis. Thus the results 
could not only be compared with those of the more restricted real MONSTER 
approach but also with those of exact shell model diagonalizations, which are not 
available in larger model spaces. It was demonstrated that the complete shell model 
spectrum of 20 Ne is reproduced exactly by the complex MONSTER approach. For 
the heavier nucleus 22 Ne still an excellent approximation to the shell model spectrum 
is obtained. So, e.g., the average deviation of the energies of yrast states from the 
exact results amounts to only 0.2% of the shell model ground state energy. In both 
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nuclei the complex MONSTER approach yields considerable improvements with 
respect to the previous real MONSTER model. As compared to the above stated 
0.2%, the real calculations give in both nuclei for the yrast energies an average 
deviation of 0.9%. 

Already the real MONSTER approach was a rather useful tool for nuclear structure 
studies. It has been successfully applied to light as well as medium heavy nuclei. So, 
e.g., a rather nice description of various nuclei in the mass 130 region was obtained 
( ||29|| , p0|). Since the complex MONSTER takes many additional correlations into 
account, we expect that this approach will develop into an even more powerful 
instrument for nuclear structure investigations in large model spaces. 

Finally, we would like to stress that even in the complex MONSTER approach pre- 
sented here, the underlying HFB transformations are still symmetry-restricted. Only 
if the additional violation of time-reversal invariance and axial symmetry is admit- 
ted, really optimal symmetry-projected vacua could be obtained by the VAMPIR 
variational procedure. In this case many additional correlations would be considered 
already on the mean field level. So, e.g., an unrestricted VAMPIR calculation in an 
lsOd-basis has 552 independent variational degrees of freedom as compared to the 56 
ones present in the complex VAMPIR mean fields. As a consequence, obviously also 
the corresponding MONSTER configuration spaces are considerably increased and 
an even more detailed description of nuclear states can be expected. Unrestricted 
calculations along these lines are planned for the future. 
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Appendix 

A Matrix elements in the "e/o" basis 

In this section it is shown that in the "even-odd" basis because of its special proper- 
ties under time-reversal only either purely real or purely imaginary matrix elements 
are obtained. As example the matrix elements of the Hamiltonian in between dif- 
ferent two-quasiparticle configurations are presented. 
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(Fa/3; SMe\H\F-f5; SMe) 
1 
2 



+ (-l) K «- K «(FaP\&_ K ^ M H& M _ I ^\F'yS) 
+ n s (-l) Is - K ^ [(FaP\e s Ka0M HQ s M _ K jF^) 

+ (-l) K °e- K ^(FaP\& Ka3M Mi K jF 1 8) 



-{(FaPie^He^jF^) 



|F 7 5) 



+ (Fa(3\Q^ M H& MKW 
+ 7rs(-l) Is - K ^ [(Fa(3\Q s KafjM m s M _ K jF^5) 

+ {Fap\e s KagM He s M _ K jF^sy 

= (Fam^HeixjF^) 

+ n s (-l) Is - K ^ 3?e (FaP\e s Ka0M HQ s M _ K jF^). 

In the same way one obtains 

(Fa/3; SMo|iJ|F 7 5; SMo) 

-7r 5 (-l)^-^ (Fa/3|e^ M iJe£^ 7 jF 7 5>, 

while for the mixed matrix elements one gets 

(Fa/3; SMe|#|F 7 <5; SMo) 
1 



(Fa/3|e^ M F,e^ 7 jF 7 5> 

_ (-l) K «f>- K -i'{FaP\i 
- ns(-l) Is - K ^[(FaP\e s KagM m s M _ K jF^)) 



_ (-l)^-^(F^|6^ flM ^_ X7 jF75) 



_ (F^|e^ a/3M M^ 7 JF 7 5) 

= i {9m (Fa/3|e| a/3MJ ff9^ 7 jF 7 5) 

-^(-1)^-^ 3m (FaP\e s KagM m s M _ K jF^5)}, 

(Fa/3; SMo|#|F 7 6; 5Me) 
= i {9m (Fa/3|e| Q ^e^ 7 jF 7 5) 

+ 7rs (-l)^-^ 9m (Fa(3\& s KafjM m s M _ K jF^)} 



and 



25 



= ? {-3m (F^S\e s K ^ M He s MK jFaP) 

= (F^5; SMe\H\Fa[3; SMo}* , (84) 



respectively. 



B Symmetries 

In this appendix some symmetry relations for matrix elements needed in complex 
MONSTER calculations are given. We restrict ourselves here to the Hamiltonian. 
For the overlap matrix elements similar expressions are valid. 

For the rotated matrix elements of type 

MF 1 ^ 2 ;^,*)]^ = (F'lHR^^xWaiFVpiF'M) (85) 

one gets 

[hxiF'F^-^-x)}^ = [h\(F 1 F 2 ; ip, x)]oap 
[h\{F 1 F 2 \ —(p, — x)] 0o/ g = —[hx(F 1 F 2 ; ip, x)}*mp-> (86) 
while for the hermitean conjugate matrix elements 

[h\(F 1 F 2 ; ip, x)]a/3o = (F 1 \a / 3(F 1 )a a (F 1 )HR('d, ip,x)\F 2 ) (87) 

one obtains 

[h\(F 1 F 2 ; —ip, — x)] Q/30 = [h\(F 1 F 2 ; ip, x)]*^ 

[h\{F 1 F 2 \ i?, —ip, —x)] a po = —\h\{F l F 2 ] <p>, x)]*apo- (88) 

If the same vacuum is used on both sides, the above two types of matrix elements 
are connected via 

[h\(F 1 F 1 ; —ip, — x)] a/ j = (-^MM^^x)]^ 
[h\(F 1 F 1 ; —ip, —x)]a0o = (—l) Ka0 [hx(F 1 F 1 \ ip, x)]oa$ 
[h\(F 1 F 1 ; —(p, —x)] a po = (-^MF'F^^x)}*^ 
[h\(F 1 F 1 ; —(p, —x)]apo = (-l)MM^ 1 ; x)W (89) 
which yields for [h 02 {F 1 F 1 \ tl)] a/3 and [hf(F l F l ; ft)] Q/3 

[h°\(F 1 F 1 ; ip, x)] a p = (-l) K « 9 K(F^,cp, X )y 
[hf{F'F^^x)U = -(-l^KiF'F^^x)}^ (90) 

26 



The same relations hold for [^(fi)]^ and [g^ 1 ^)]^. 
For the rotated matrix elements 

[h\(F 1 F 2 ; i?, </?, x)]afSj5 = 

(FV^F^a^F 1 )^^^^)^^ 2 )^!^ 2 )!^) (91) 

we get 





-<P, 


afiyS = 


[hx{F^F 2 -^ X % m 




-<P, 




[h\(F 1 F 2 ; ip, x)]*ap~fS 


[h\(F l F 2 ] 


-V, 


~~ XjlaPfS = 


— [h\(F 1 F 2 ; ip, x)]*apys 


[hxiF'F 2 ;^, 


-<P, 


_ ^)]a/37<5 = 


~[h\(F 1 F 2 ; ip, x)]*apys 


[hxiF'F 2 ;^, 


-<P, 


~~ X)] a f3-y8 = 


~[h\(F 1 F 2 ; ■&, ip, x)]*a/3y5 


[hxiF'F 2 ;^, 


-<P, 


~ X)] a f3jS = 


~[h\(F 1 F 2 ; ip, x)]*ap*f5 


[hxiF'F 2 ;^, 


-<P, 


~ X)] a 0^S = 


[h\(F 1 F 2 ; <p, x)]*a/3js 


[hxiF'F 2 ;^, 


-<P, 


~~ X)] a/ 3 7 <5 = 


[h\(F 1 F 2 ; i?, ip, x)]*a/3j8- 



For identical vacua on both sides furthermore 

[h\(F 1 F 1 ; — </?, —x)] a p-yS = {—l) Ka0+K ~ lS [h\{F 1 F 1 \'&, (p, x)]* 5q/3 (93) 

holds. 

Finally, for the Hamiltonian matrix elements in between two one-quasiparticle con- 
figurations 

\h x (F l F 2 -$,^x)U = (F'MF^HR^^^W^F 2 )^ 2 ) (94) 

the relations 

[hx{F 1 F 2 \ i9, —p, —x)] a p = [hx{F l F 2 -^^x)U 

[hx{F 1 F 2 \ —ip, —x)] a p = —[h\(F 1 F 2 ; ip, x)]*ap (95) 

are found. 
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Figure Captions 



Fig. 1 The energy spectrum for the nucleus 20 Ne as obtained by three different 
approaches : the complete shell model diagonalization (SCM, solid lines), 
the complex MONSTER (CM, dashed lines) and the real MONSTER (RM, 
dotted lines). The energy is given relative to the 16 core. Only excitation 
energies up to f 6 MeV above the shell model ground state are presented. The 
real MONSTER calculation has been based on the corresponding + VAMPIR 
transformation. The complex MONSTER results, which do reproduce all shell 
model states exactly, are independent of the particular choice of the underlying 
transformation. 

Fig. 2 Same as Fig. 1, but for 22 Ne. For each spin only the lowest shell model 
states are presented. The complex MONSTER was here built on the complex 
VAMPIR solution obtained for the + ground state. 

Fig. 3 The energies of the yrast states in 22 Ne as obtained by various approaches : 
the shell model (SCM), the complex MONSTER (CM) (again on top of the 0+ 
VAMPIR solution), the complex VAMPIR (CV), the real MONSTER (RM) 
(based on the corresponding + VAMPIR solution, too), and the real VAMPIR 
(RV). In the last approach only the even spin states are accessible. Spin and 
parity are indicated on the l.h.s. of each level. Again the energy is given 
relative to the 16 core. 

Fig. 4 The five lowest 0+ MONSTER states of 22 Ne are plotted versus the spin of 
the VAMPIR transformation which was used in each calculation. Solid lines 
refer to the complex MONSTER, dotted lines to the real MONSTER results. 
In the latter case only transformations with even spin values could be used. 
For comparison the lowest three + shell model (SCM) energies are also given. 
Again the energy is given relative to the 16 O core. 
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